CERN-TH. 741 1/94 
NTUA. 10/94 
IOA.307/94 



LOW ENERGY THRESHOLDS AND THE SCALAR MASS 
SPECTRUM IN MINIMAL SUPERSYMMETRY 



G. K. Leontaris* 

CERN 
Theory Division 
1211 Geneva 23 
Switzerland 
and 

N.D. Tracas 

Physics Department 
National Technical University 
GR-157 80 Zografou, Athens, Greece 

ABSTRACT 

We discuss low energy threshold effects and calculate the sparticle masses 
in the context of the Minimal Supersymmetric Standard Model. We pay 
particular attention to the top squark and the Higgs mass parameters, and 
calculate the top Yukawa corrections, taking into account the successive de- 
coupling of each particle at its threshold. We discuss the phenomenological 
implications in the context of the radiative symmetry breaking scenario. 
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The Minimal Supersymmetric Standard Model [] (MSSM) has by now 
been accepted as the most natural extension of the Standard Model of Strong 
and Electroweak interactions. As the recent experiments approach closer 
to the energies where some of the superpartners seems to acquire their masses, 
it is very important to obtain higher precision in the theoretical predictions 
of the scalar mass parameters, Yukawa coupling corrections, threshold effects 
etc. In recent analyses j|, ||, |^, 0, |[, it has been shown that the semi-analytic 
procedure in the calculations of the above quantities can offer the possibility 
of investigating reliably the above effects. Moreover, the advantage of ana- 
lytic expressions for the low energy parameters is more than obvious: one can 
extract easily information about the role of the input parameters at the GUT 
scale (mo, my2, h ta ), since in the analytic procedure the low energy mea- 
surable quantities can be expressed in terms of calculable functions of the 
former with the boundary conditions incorporated into these expressions. 
Nevertheless, detailed theoretical predictions may be still pushed further by 
estimating higher order effects, threshold corrections, etc. 

In recent works it has been shown that the unification scenario |1 survives 
even when various uncertainties arising from several sources (GUT and low 
energy SUSY thresholds, input values for coupling constants, experimental 
uncertainties etc.,) are taken into account^, || ||, 0- For example, in ref. [§]], 
it has been shown that an effective low energy SUSY scale can be defined 
which, for a realistic mass spectrum can account for low energy threshold 
effects. More recently a more accurate way to estimate the uncertainties 
of such effects, which may also take into account two loop corrections, has 
been explored. 

Predictions of low energy parameters turn out to be very sensitive in all 
the above mentioned threshold effects. In estimating these effects, it has 



been shown fl2|1 that it is adequate to use the "step" approximation in the 
definition of the beta function coefficients. Thus the values of the weak mix- 
ing angle (sin 2 9w), the electromagnetic (a em ) and the strong (03) couplings, 
and other low energy parameters can be given by their one loop formula 
with the addition of two small correction terms arising from two-loop and 
threshold effects in a semi-analytic procedure^. In calculating however the 
scalar masses themselves, one should be careful in particular for those affected 
from the top- Yukawa coupling h t . In this case due to the non-negligible con- 
tribution of h t , the evolution of m~ squark and m# 2 mass parameters are 
determined by a coupled differential equation system. Thus in addition to 
the successive changes of the gauge coefficients as the sparticles decouple, the 
boundary conditions at each sparticle's threshold should be treated carefully. 

In the present analysis, we wish to investigate the above effects in the 



^or reviews see for example [fij 
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context of the MSSM assuming that the gauge couplings unify in a simple 
non-abelian gauge group at an energy scale close to 10 16 GeV. We find it useful 
to adopt a semi analytic procedure and provide specific formulae for all the 
involved parameters, and compare them with those of previous estimates 
where such effects were not taken into account. 

In particular the following issues will be discussed. We will start assum- 
ing that the radiative symmetry breaking (RSB) scenario [0 is an effective 
mechanism operating in the usual sense at low energy, i.e. by driving one of 
the Higgs mass-squared parameters negative at energies close to mz- We are 
going to use one-loop corrections to the effective potential and estimate the 
effects in the -parameter which plays an essential role. Next we calculate 
the exact contributions of the trilinear parameter A and compare the results 
with previous estimates where these corrections were not included. Finally, 
we are going to calculate the scalar masses for various choices of the initial 
values mo,TOi/2 taking into account the afore-mentioned threshold effects. 

Starting at the GUT scale with a particular gauge group, one chooses 
specific values for five independent parameters, namely m , mi/2, /•*, A and 
B. In the simplest case, all the scalars have a universal Q mass m$. The 
masses evolve down to low energies where one expects that one Higgs mass- 
squared parameter becomes negative. This triggers the SU(2) x U(l) sym- 
metry breaking. The calculation of the mass-squared parameters needed to 
check if this scenario is valid, requires the solution of the coupled differential 
R.G. equations obeyed by these parameters. 

In the case of the small tan f3 ~0(1) scenario (tan /3 is the ratio of the two 
v.e.v's), one may approximate the relevant differential equations as follows 
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where only the top quark Yukawa coupling h t has been kept. The coefficients 

2 for recent discussions where non-universal conditions at the GUT scale are assumed 
see 111 
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q are given by 



while Mj's are the gaugino masses and t = \nQ. 

The differential equation for Hi can be solved straightforwardly, since 
it is independent from the others. The remaining four differential equations 
define a coupled system which depends strongly on the top Yukawa coupling. 
Making the identifications 

mi L = rhi , mi R = m 2 and m H , 2 = m 3 

the solution of the system is found to be 



m 2 Hl (t) = m 2 + C 3 (t)mi /2 (6) 
m 2 n {t) = ml + C n {t)m\ /2 - n5 2 m (t) - n5 2 A (t) (7) 
A(t) = q-\t)(A G -mi /2 I A {t)) (8) 

where all the relevant functions are presented in the Appendix. 

The scale dependence of the Yukawa coupling h t (t) can be found by solv- 
ing the RGE for that coupling 

dh t 1 / ,o A o\ , , .13 „ 16. 



dt 8n 2 ^-gc^j^, where Q = {-,3, y } (9) 

with the well known solution, 

htit) = h tGlu (t)q- 1/2 (t) (10) 

In all previous equations, the subscript G denotes the corresponding value 
at the GUT scale. 

There are four arbitrary parameters entering the above formulae, namely 
n^o, m i/2, A G and h tG . The first three of them, as was already pointed out, are 
the soft input mass parameters at the unification scale M G . Since the above 
solutions enter the minimization of the Higgs potential V(Hi, H 2 ), their range 
can be phenomenologically constrained by the requirement of generating a 
stable minimum for this potential. Of course, a crucial role is played also 
by the fourth parameter, h ta , which should be large enough to drive the 
i/2-Higgs mass-squared parameter negative and give a phenomenologically 
acceptable vacuum. 
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Experimental evidence |HJ as well as theoretical expectations ||16|| treating 
the Yukawa couplings as dynamical variables, indicate that the top mass 
requires a large top- Yukawa coupling, close to its infrared fixed point, i.e. 
(m t / sin (3) ~ 190 GeV. Therefore, since only the ratio m t j sin (3 enters in the 
relevant running scalar mass parameters, one may safely conclude that in 
the fixed point solution for the top-mass, their values depend mainly on the 
initial values mo and mi/2- The above argument may be more transparent 
if one writes the above mass-formulae in the limit of the infrared fixed point 
of h t . One then obtains M 
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The contribution of the trilinear parameter A introduces one more input 
parameter at the GUT scale, but its role is less significant as long as Aq is 
of the order of m (as expected). Indeed, by writing Aq = A mo, after some 
algebraic manipulations one can show that S A is written as follows 

= ^i>' )d K) = -w)L" V)t W- w (12) 

Expanding A 2 (t), using Eq.(H), we write formally 

5 2 A = 5 2 Ai A 2 ml + 5 2 A2 A m m 1/2 + 5 2 Aa m 2 1/2 (13) 



where 



^ = -IMw)- 1 ) (14) 



6q(t) \q(ty 



= \— Ut)-lA(t)\ (15) 



= (^) J ^ - 2 ^(*)M*) + 2/^)) (16) 



Ut) = f q(t')C A (t') lA (t')dt> (17) 



7x(t) = W = Er — -i (is) 



The coefficients 5a, depend on simple integrals of scale dependent pa- 
rameters. Evaluation of the relevant integrals give the following results (for 
m t ~ 175 GeV and SUSY breaking M s = (500 - 1500) GeV) 

S\ = 0.005^mo - 0. 020 A Q m Q m 1/2 + (0.146 - 0.155)m? /2 (19) 
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which add too small corrections to the solutions where A was ignored. More- 
over, these corrections become even smaller |5], [7| as m t approaches its infrared 
fixed point value. 



In the above equation, the simplified assumption was made that all scalar 
masses decouple at the same scale, namely Ms- In the running of the RGE's it 
is assumed that there is a great "desert" between the GUT and the weak scale 
while a low energy SUSY scale Ms is assumed so that for energies lower than 
Ms one uses the Standard Model beta function- and c, coefficients. However, 
as has been already pointed out, a more careful treatment should also take 
into account threshold effects due the successive decoupling of these scalar 
masses from the spectrum at different scales. In the semi-analytic approach 
of ref. [§J an effective scale Mg was assumed which can be roughly estimated 
to be 



to account for the SUSY scalar mass effects. The parameter can also be 
given in terms of known parameters by solving the minimization conditions of 
the neutral Higgs potential. Taking also into account one-loop contributions 
to the superpotential, we may obtain the following approximate formula for 



the tree level contribution while 7] and Q are defined in the Appendix. For 
tanfi > 1.1, is less than 1.5 GeV. One therefore, for sensible values of the 
parameter \fi\ < (1 — 2) TeV, could define a reasonable scale Mg ; below that 
scale the beta function coefficients should turn to their non-supersymmetric 
form. 

The precise effects, however, are found by the successive change of all the 
beta function dependent coefficients at each particle's threshold. Assuming 
only one-loop corrections, in the Minimal Supersymmetry with three families 
and two Higgses, one can write the b[s in the following form[[|, [L7], [Tj| 




(20) 
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In the above formulae for b^s, H stands for the higgsino contribution W 
for the winos, etc., while for any particle's threshold with mass m 2 ., we have 
denoted 9 Si = 6(Q 2 - w 2 .). 

In our semi-analytic approach, when evolving the gauge and Yukawa cou- 
plings as well as the scalar mass parameters down to low energies, we find 
it sufficient to define the following ^-changing scales: We assume a common 
scale Ql for the decoupling of UL 12 ,dL 12 sparticles, while we assume that 
they are not very different from their mass eigenstates. The next scale is the 
one defined as an average scale Qr of their right-handed partners. In the case 
of the universal scalar masses at the GUT scale, however, these two scales 
could not differ substantially unless mi/2 is very large (> ITeV). Therefore, 
the thresholds arising between these two scales are not expected to have a 
significant effect. 

We define as a third scale (subsequently denoted with t\ = Infhi) where 
we change the bi and q coefficients, the scale where the t^-squark acquires 
its mass. As it is expected, due to the large negative contributions from the 
top-quark Yukawa coupling, this mass should be substantially smaller than 
those of the L, R squarks of the first two generations. Finally, we define 
two more new scales above the weak and top mass scale, namely the scale 
were the tR squark gets its mass (subsequently denoted with t 2 = Inrh^) and 
an average scale for all other contributions (sleptons etc). The hierarchy of 
these two latter scales depends strongly on the point of the parameters space 
(m ,m 1 /2) one has chosen. A simple inspection of the obtained evolution 
equations for their masses shows that slepton masses are larger than mz for 
m values much bigger than mi/2 while the opposite is true when mi/ 2 > m Q . 
Of course all the above scales should be carefully incorporated in the analytic 
formulae presented above. All the relevant integrals (see Appendix) should 
split into sums over the various scales. 

Our next step is the determination of the mz L and mx R mass parameters 
which also define the afore-mentioned scales where the beta function coef- 
ficients should also change. For a given (mo,mi/2) pair, one can compute 
the left and right squark masses of the first two generations. The negative 
corrections <5^(ti) may then determine the mass of the ti squark. Below this 
scale, top- Yukawa negative corrections should not include contributions from 
diagrams involving ti- Therefore in the range defined by mz L > Q > mz R (i.e. 
for the range (£i,£ 2 )) the evolution equations for m 2 = rh\ and m 2 H2 = m| 
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can be written as follows 

rh 2 2 = Tr^{t 1 ) + C 2 (ti } t)ml /2 -2S{t u t) (25) 
m 2 3 = rr^(ti)+C 3 (ti,t)m? /2 -3<y(ti,t) (26) 

where m^(ti) are the mass parameters calculated at the scale t\ = ln{m^ L ). 
Top- Yukawa corrections 8(ti,t) contain now the sum only of tn squark and 
the H 2 Higgs mass parameters, and possibly (depending on the specific values 
of the (mo, 7711/2) pair) the A(t) trilinear parameter, 

S(h,t) = £^r (^ml{t>)+0(t> -t A )A\t')}dt' (27) 

where defines the logarithm of the scale at which the trilinear mass pa- 
rameter stops running. The above corrections can be calculated easily, by 
solving (|25|), fl26|). The result is 

S(h,t) = -q-Hh,t) [%l( tl ,t')u (t')^dt' (28) 

Jti 07T Z 

where 

uo(t) = 2m 2 + [C 2 (ti, t) + C 3 (ti, t)] m\ /2 - 5[5 2 m (h) - 5 2 A (h)] (29) 

After the decoupling from the massless spectrum of the squark at the 
scale t 2 = ln(mi ), one ends up with only the Higgs mass parameter whose 
evolution from t 2 until the transmutation scale is given by the formula 

ml(t) =ml(t 2 )+C 3 (t 2 ,t)mj /2 -3q-^(t 2 ,t) j\l(t 2 ,t')v Q {t') } ^-dt' (30) 
where 

v Q (t) =ml(t 2 ) + C 3 (t 2 ,t)m 2 1/2 (31) 

In tables / and // we present the calculated scalar mass spectrum for two 
initial values of the h ta , (case //, very close to its infrared fixed point value), 
and representative choices of mo, mi/2 pairs. The last column of these tables 
shows the prediction for a s low energy parameter, for each sparticle spectrum. 
In obtaining our results, we have worked in the low tan/3 regime while we 
have allowed 10% deviations from the GUT relation h^ta) = h T (t G ). On 
the other hand, the obtained values for sin 2 Qw and m t are consistent with 
the relation sin 2 6 w (m z ) = 0.2324 - 10~ 7 x {(m t /GeV) 2 - 143 2 } ± 0.0003. 

As it can be inferred from the tables, the ti,R squarks and the average 
sparticle spectrum is lighter for larger top couplings and top mass. When 
mo is relatively small, slepton masses are the lighter sparticles while mi R 
becomes light when m\i 2 mo- Of cource, right sleptons (not shown in the 
tables) have slightly smaller masses than their left partners. 
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Table I. The supersymmetric sparticle spectrum together with the 
corresponding a s prediction in the low tan/3 scenario (tan/3 ~ 1.5) and 
m t < 165 GeV, fixing sin 2 9w around its central value ~ .232 

Now let us discuss the effects of the decoupling of the havier quarks in the rest 
of the sparticle spectrum. The successive decoupling of each scalar mass term 
from the relevant differential equation, has modified the negative corrections 
induced by the top-coupling below the ti squark mass. In the low tan(3 
scenario this treatment has a direct effect only on the t^-squark and the m^ 2 
mass parameters. In particular, in our semi-analytic treatment we observe 
that the £r mass has increased by (1 — 5)% , relative to a naive treatment 
of the boundary conditions at each particle's threshold 0, depending on the 
specific choice of the (mo,mi/2) point. Such corrections are therefore of the 
same order but with the opposite sign of the A-trilinear term corrections 
given by the formula ([S]). The maximal change occurs when the top Yukawa 
coupling gets closer to its fixed point value {Table IT). The effect however 
is still small, since the modified equation fl2"5|) runs only on a small region, 
namely between the m th and m tR scales. A larger effect is found in the m,H 2 
running mass whenever the transmutation scale is substantially different from 
m tL . 

Finally, as it has been pointed out above, the t^ mass is smaller than 
the slepton masses in specific regions of the (mojin^) parameter space. A 
qualitative picture of the mj and mj variation in terms of tuq and is 
given in the Figure. 
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Table II. The supersymmetric sparticle spectrum together with the 
corresponding a s prediction for h tG Yukawa coupling close to its 
non-perturbative value, and m t ~ 170 GeV. Again sin 2 9w has been fixed to 

its central value ~ .232 

We can see that there is a considerable fraction of the parameter space 
(mo, mi/2) which allows solutions of relatively small mi R . For example, in 
the last entry of table II, mi is of the order of the top-quark mass. This 
would imply that, after the diagonalisation of the squark mass matrix, the 
light physical mass eigenstate m^ could be as small as 150 GeV. This gives 
hope that future experiments may discover supersymmetric signatures. 

To summarize our results, we have used a semi-analytic approach to cal- 
culate the supersymmetric spectrum in the small tan/5 regime, taking into 
account low energy threshold effects. We have given special emphasis to 
top squark and Higgs mass parameter calculation, which in the presence of 
a heavy top quark receive large negative contributions. We have examined 
in detail the effects of the 'decoupling' of the heavier sparticles from the 
renormalization group equations of the lighter ones, and we have found that 
our treatment of the boundary conditions, results in an increase (1 — 5)% of 
the tR mass parameter, compared to a naive treatment. Corrections on the 
scalar masses from the trilinear parameter A are treated also analyticaly and 
found to be of the same order for moderate initial values (Ac ~ — V3m ) 
Furthermore, we have examined general properties of the sparticle spectrum 
and observed interesting correlations. Thus, large values of m compared to 
that of mi/2 imply that mi R is lighter than the left slepton masses, while the 
opposite is true for > m . Moreover, in the small tanfl regime that we 
are examining here, for a considerable fraction of the (m , 7771/2) space, a light 
t-squark (~ 150 GeV) can be obtained, which might be found in accessible 
energies by experiment in the near future. 
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APPENDIX 

The scale dependent coefficients in the scalar mass solutions (|6|,[7||f) are 
given by the following general formula 

3 n 

C n(h, t) = £ —V - ocKti)) , n = 1, 2, 3 (32) 

1=1 Z °i a i a 

with the identifications C\ = Cq, C% = Co and C3 = Ch 2 - We define 
C n {t) = C n {t G ,t). 

The gauge dependent functions ju(t), and g(t) which arise from the top- 
Yukawa differential equation have the form 



^)=nn(#r] (33) 



n i= 



3 T / " 

3// 7. . . 3/iL /* 

where the index n runs over all the intermediate scales. Again we define 
q(t)=q(t G ,t). 

The negative Yukawa contributions of Eq.(^) are found to be || 



«( ( 1 . ( ) = l + ^f^) = l- 1 #/^(0*' (34) 



«,(*) 



(2£^)W« + ™WM). (35) 

where v = 246 GeV, while the quantities I, J, I a are functions of scale de- 
pendent integrals given by: 



m = - 1 lUt'w, j(t) = - e / lUtmt'w (37) 

A M = f G ^A\t')dt\ E\{t) = (t')Ai(t')^' (38) 

Jt 07T Z Jt 

I A (t) = f q{t')C A {t')dt> = ^jZcfa lG f q (t')^P-dt' (39) 



The minimization conditions of the tree-level neutral Higgs potential give 
the following solution for the |/io|- parameter M 



1 rr + 2 , / k 2 J \ 9 9 1 V2 

{fc^T m ° + (p^T7 " l ) m ^ ~ M| ) (40) 



V2 
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with k = tan (5. 

The parameters 77, Q entering the one loop formula are given by 
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with p = m t /Mz, p =< mi > /Mz and po the tree level parameter defined 
in ([H]). Finally the t-squark mass-combinations M\ L ± M\ R are given by 



M\ L + M 2 
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